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ABSTRACT 

We discuss the galaxy-galaxy-mass three-point correlation function and show how to 
measure it with weak gravitational leasing. The method entails choosing a large of 
pairs of foreground lens galaxies and constructing a mean shear map with respect 
to their axis, by averaging the ellipticities of background source galaxies. An average 
mass map can be reconstructed from this shear map and this will represent the aver¬ 
age mass distribution around pairs of galaxies. We show how this mass map is related 
to the projected galaxy-galaxy-mass three-point correlation function. Using a large 
N-body dark matter simulation populated with galaxies using the Halo Occupation 
Distribution (HOD) bias prescription, we compute these correlation functions, mass 
maps, and shear maps. The resultant mass maps are distinctly bimodal, tracing the 
galaxy centers and remaining anisotropic up to scales much larger than the galaxy sep¬ 
aration. At larger scales, the shear is approximately tangential about the center of the 
pair but with small azimuthal variation in amplitude. We estimate the signal-to-noise 
ratio of the reconstructed mass maps for a survey of similar depth to the Sloan Digital 
Sky Survey (SDSS) and conclude that the galaxy-galaxy-mass three-point function 
should be measurable with the current SDSS weak lensing data. Measurements of this 
three-point function, along with galaxy-galaxy lensing and galaxy auto-correlation 
functions, will provide new constraints on galaxy bias models. The anisotropic shear 
profile around close pairs of galaxies is a prediction of cold dark matter models and 
may be difficult to reconcile with alternative theories of gravity without dark matter. 

Key words: gravitational lensing 


1 INTRODUCTION 


The last few years have seen the emergence of an observationally consistent cosmological model based on firm physical ideas. 
The hot big bang model of the expanding universe, whose structure and dynamics can be described by general relativity (GR), 
has a number of cosmological parameters that are now well constrained by n umerous observations. WM AP measurements 
of the cosmic microwa ve background radiation (CMB) JSnergel et a i] 12001^1. supernova measurements jRiess et al.l llQQ^ : 
IPerlmutter et IDB , and measurements of large scale structure et among others, have settled on a 

concordance model of a spatiall y flat universe with m atter density about 25% of critical . This concordance mo del is consistent 
with the ages of the oldest stars il.Timenez et alll99fitl , Hubble paramet er measurements (IPreed man et al ■1200 ill . and primordial 
abundances of light elements from Big Bang Nucleosynthesis (BBN) llBurles et alJl200lll . 

One of the main requirements of a successful cosmological model is that it connect measurements of the early universe, 
e.g. the CMB and BBN, to the measurements of the local universe of galaxies and stars. The large scale structure of galaxies 
that we observe in galaxy surveys must have formed through gravitational collapse of the small fluctuations left over from 
that earlier time. The properties of these large-scale mass densities can be predicted from the initial conditions as observed in 
the CMB combined with our understanding of gravity as described by general relativity. However, since the Universe in our 
model is dominated by dark matter, we need to make the connection between the distribution of galaxies relative to the dark 
matter if we want to compare these observed structures to what is predicted. 

The relation between the galaxy density and the dark matter density is known as the bias. For the concordance model to 
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be in agreement with the data on galaxy clustering, it is required that on large scales the bias parameter (which we will define 
shortly) is close to one. On smaller scales the bias must be different from one but of order unity. Moreover, the relative bias of 
different galaxy types is known to depend on spectral type and luminosity. One can do better than simply treating the bias as 
a free parameter. Galaxy biasing is after all closely linked to the physics of galaxy formation. Although galaxy formation is not 
completely understood, it is becoming increasingly amenable to theo retical modeling and numerical simulation. G as dynamical 
N-body simulations and Semi-analytic models of galaxy formation llRensnu et alJl^nnl: iKanffmanu et each predict 

how galaxies are biased with respect to mass. Nevertheless, sufficient uncertainty on the details of galaxy formation remain 
that observational constraints on the bias are needed. 

The apparent complexity of the bias has motivated the development of new probes of large-scale structure and bias that 
do not rely on luminous galaxies as tracers of the underlying mass. Gravitational lensing is a prime example. Tensing is ideally 
suited to studying dark matter since it is directly sensitive to the total mass density along the line of sight. Other probes of 
dark matter such as dynamical measurements must rely on some assumptions such as that a system has reached dynamical 
equilibrium. By contrast, lensing provides a clean probe of the dark matter. 

Galaxy-galaxy lensing (GGL), in particular, probes the dark matter density profile around luminous galaxies. In galaxy- 
galaxy lensing, one measures the mean tange ntial distortion of a large number of source galaxies behind foreground lens 
galaxies by stacking the signal from the latter llSheldon et al.lEo04ll . To date, lensing measurements and theory have focused 
on the two-point galaxy-mass correlation function (2PCF). This correlation function can be measured from galaxy photometric 
surveys and compared to the results of N-body dark matter simulations. These measurements provide direct constraints on 
bias models. Although they do not measure the galaxy bias parameter, b , directly, they do measure a closely related quantity 
and so these measurements do test the bias model. 

While the 2PGF is the lowest order measurement of how much a distribution is clustered, it does not carry complete 
statistical information about a non-Gaussian distribution. The rest of the information is encoded in the infinite set of N-point 
correlation functions. In particular the higher order (N > 2) functions provide the information about the coherence and shapes 
of large-scale structures such as filaments, walls, and voids that are apparent in both data and simulations. These higher order 
correlation functions add further constraints to be met by cold dark matter (CDM) models and new probes of the bias. 

For the most part observational studies of higher order correlation functions have measured correlations between galaxies. 
In this paper we will describe a method for measuring the galaxy-galaxy-mass three-point correlation function (GGM3PGF) 
with weak lensing. This method is a generalization of the method for measuring the galaxy-mass two-point correlation 
function (GM2PGF) with galaxy-galaxy lensing. Just as the GM2PGF provided new constraints on the bias model, so too 
does the GGM3PCF. While the GM2PCF is a measure of the isotropic mass density around galaxies, the GGM3PGF gives us 
information about the shape of the mass distribution around pairs of galaxies. Together with galaxy N-point functions and the 
GM2PCF, the GGM3PGF provides independent constraints on models of galaxy biasing that are crucial for understanding 
both large scale structure and the processes of galaxy formation. 

It is reasonable to assume that GDM models should predict that the mass around a pair of galaxies will have an anisotropic 
distribution. A priori it may be purely elliptical or clumped tightly about each galaxy with a distinctly bimodal appearance. 
This expectation comes from the fact that halos are not themeselves spherically symmetic and that the galaxies should also 
trace this ellipticity. Also the galaxies themselevs have mass associates with them; stellar mass at the least and probably 
dark matter sub-halos. Whatever the details of the distribution is, it will create a gravitational lensing effect that is also 
anisotropic. This anisotropic lensing signal should be measurable with current and future weak lensing surveys. 

Another use of the G GM3PGF is to test for the existence of dark matter itself. Alternative theories of gravity with no 
dark matter, e.g. MOND llMilgromlll QS.jh . must eventually make prediction for gravitational lensing. This predicted lensing 
profile must agree with the lensing observations. With no dark matter, the lensing effect is only caused by the visible matter. 
Two galaxies that are close together will simply look like a monopole at large distance from them and the lensing effect 
produced should be nearly isotropic regardless of the form of the lensing profile. However, dark matter models predict (as we 
will demonstrate) that galaxy pairs have an extended dark matter halo about them that is anisotropic even at large distance 
from their center. This causes an anisotropic lensing effect at these large distances. Thus, dark matter models and alternative 
models should make different predictions for the lensing anisotropy and thus test the idea of the existence of dark matter. 
Hoekstra et al. 2003 have already claimed to have measured halo flattening around single galaxies, a similar effect, with 
weak lensing and argue that it would be difficult to reconcile this anisotropy with theories such as MOND. Measurements 
of the GGM3PCF would allow one to extend this result to scales larger than the size of individual galaxies. Also, since the 
GGM3PCF deals with only the locations of galaxies, it would not have to rely on simulations that predict how the shape of 
the galaxy light correlates with the shape of the halo. 

In this paper, we show how the lensing pattern around pairs of galaxies is related to the GGM3PGF. In particular, 
we show how the 2D projection of the 3D GGM3PGF is directly measur able with weak lensing measurements such as 
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the measurements being conducte d 


___ iMc.Kav et alJl^Olli with data from the Sloan Digital Sky Survey 

llVork fc SPSS collaborationlboOfll . We use an N-body dark matter simulation populated with galaxies using a bias model 
chosen to reproduce the results of the hydrodynamical simulations and make numerical calculations of the various two and 
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three-point correlation functions in both 2D and 3D. We calculate the lensing effect and estimate the amount of noise that 
should be present to study the signal-to-noise that can be obtained. We conclude that the GGM3PCF can now be measured 
with the present SDSS data. 


2 CORRELATION FUNCTIONS AND THEIR APPLICATIONS TO COSMOLOGY 


We begin by reviewing galaxy correlation functions and then describe their extension to galaxy-mass cross-correlations. The 
probability of finding a galaxy (or any mass point) in an infinitesimal volume of space dV\ is simply dPg = PgdVi where pg 
is the mean galaxy density. The probability of finding one galaxy in dVi and another in dVi defines the two-point correlation 
function (,gg{r) 

dPgg = pldV^dV2 [l + ?99(r)l 

The function CgaG) can be thought of as the contribution of clustering to the joint probability in excess of random. By the 
assumed large scale homogeneity and isotropy, ^gg can depend only on r and not on angle or location in space. The marginal 
probability of finding a galaxy in dV 2 given that one already has found one in dV\ is dP{g 2 \gi) = p^dVidV 2 + Pg^Vi = 

pgdV2 [1 + Css(?■)]• 

We can further consider the probability of finding a third galaxy in infinitesimal volume dVs in addition to the two in 
dVi and dV 2 . The resulting joint probability defines the three-point correlation function (3PCF), Csss i 

d-Pggg ~ pgdV \ dV 2 dV 3 [ 1 -f Css(^12) + Css(’'13) + Css(’'23) + Csss] 


Here rij refers to the separation between galaxies i and j. The point of defining Csss ®C) as to result in a 5-term expression 
is similar to the reason we defined a 2-term expression for the 2-point case. If the third galaxy were simply put down at random 
or not clustered with the first two, or very far from the first two, we would have dPggg = -I-Css(’' 12 )]■ Thus the 

other three terms. Css (’' 13)1 Css (’' 23 ), and Csssi only contribute the excess over this random occurrence and by symmetry of 
this argument we also get the next two terms Css (’’is) ^cid Css (’' 23 ) • The last term Csss needed since there is no reason to 
assume that the first four terms completely describe the joint probability. 

The volume averages of the 2PGF, 3PCF, and higher order correlation functions are the moments of the probability 
distribution of the smoothed density field, e.g. the variance, skewness, kurtosis etc. Since a Gaussian distribution is charac¬ 
terized just by its 2PGF, the 3PCF is the lowest-order N-point statistic that characterizes non-Gaussianity of the density 
field. Generally non-linear gravitational effects generate a non-Gaussian density distribution even if the initial density field is 
Gaussian. 

One can write down expressions for the joint probabilities of still more galaxies and these in turn will define the N-point 
correlation functions. The 3PCF Csss) due to the large scale homogeneity and isotropy, is only a function of the three triangle 
sides ri 2 ,ri 3 , r 23 (or any other parametrization of the triangle). 

One can also define cross correlation functions between different species of objects, for example, between galaxies and 
dark matter mass particles. We can write down joint probabilities of finding a galaxy and a mass particle, a galaxy and two 
mass particles, and two galaxies and a mass particle respectively (taking all the dV to be equal in size for simplicity). 


HP — 

gm — 

fiP — 

gmm — 

HP — 

ggm — 


PgPmdV 

PgPm^^ 

plpmdV^ 


[1 + igm.{r)] 

[1 -f ^ gm . Gg ^ rni '} -f ^gm (^ 9 , 7 ^ 2 ) ^ 7717 ^(^ 771 ;^, m2 ) Cffmm 

[1 + Cggirgi , g2 ) + Cswi (?' gi , m ) + ^ gm { rg2 , m ) -f Cggm ] 


( 1 ) 

( 2 ) 

( 3 ) 


These joint probabilities define the three two-point functions, ^gg , ^gm, Cmm and the four three-point functions (^ggg , 
Cggm, Cgmm and 

One can also make the same definitions in 2D for a set of points that have been projected over the third dimension. In 
this case we will refer to these 2D correlation functions as, for example, Wgm for the two-point and Zggm for the three-point 
functions. 

The uses of correlation functions in cosmology are many. In the inflation scenario, structure formation is seeded by 
random quantum fluctuations. In most models these initial fluctuations are predicted to be Gaussian distributed and therefore 
completely described by their two-point correlation ^(r) or its Fourier transform the power spectrum P{k). In this scenario, 
the 3PCF and higher order correlation functions that arise are entirely due to the onset of non-linear gravitational collapse and 
so should be predictable with dark matter simulations given a set of co smological parameters and a galaxy bias prescription. 


law, in agreement with predictions that use N-body simulation combined with HOD bias models (see Section]^. 


The first measurement s of the galaxy two-point correlation function IITotsiiii fo Kiharall QfilltfPeehle.Jl 974al : lGroth fc Peebles! 

rneull 9731 found that it was well describ ed by a power law, ^gg(r) = (r/ro) with 7 = 1.8 and ro « 5/i” 


Mpc. The most recent measurements JZehavi et al ] l2003ll are now sensitive enough to detect small deviations from a power 
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The first measurement of the 3PCF for galaxies was bv IPeebles fc GrothI For these early correlation function 

measurements, there was no redshift information, so these were measured in terms of angular coordinates 612 rather than 
physical length scales ri2. Constraints on the 3D 2PCF and 3PCF came from inverting the projections. From this first 
measurement, it was found that to experimental accuracy the 3PCF could be written in a reduced form relating it to the 
2PCFs (on scales less than a few Mpcs) as 

■*999 ~ Q + ^r)gg(0i2)u)gg(023) + U)gg(0i3)u)gg(023)l (4) 

with Q being a constant « 1. This scaling relation is referred to as hierarchical scaling. Although much effort has gone into 
deriving such an equation from first principles, no one has shown analytically why this form should be true. However N-body 
simulations do indicate that the 3PCF approximately takes this form in the non-linear regime probed by these early measure¬ 
ments. On larger scales (r > lQh~^ Mpcs) where one can treat gravity perturbati vely, this factor Q is indeed independent 
of scale however there is a non-vanishing dependence on triangle shape llFrvlIlQS^l . This large scale shape dependence is a 
quantitative statement that co-linear triangles or filamentary structures are more likely than isosceles triangles or isotropic 
structures and physically relects t he fact that gravitational collapse happens anisotro pically with velocity flows along density 
gradients (iRernardea.u et ^ llsooil . A visual inspect ion of recent N-body simulations l lJenkiiis et aljITQQsl: lEyrard et all boodi 
as well as real data from galaxy redshift surveys liHuchra, et 31.11799131 : IPeacoc.k et al .I bnnil: iTee-mark et a.lJbnO.lh show the 
existence of large-scale filaments on the outskirts of giants voids. 

The basic hierarchical scaling of Equation 2] was confirmed bv iFrv fc Seldneil (Il982ll however these early measurements 
from the Lick galaxy catalog did n ot extend to large enough s cales to reveal the exp e cted c onfiguration-dependent signa- 
ture from gravitat i onal in stability llFrvlll 98^ . Later results hy jFrienumfc (T^afiagal ill 99!Jl using the APM survey data, 
IScoccimarro et all i2001ah using the IRAS 1.2 Jy survey a.nd lFeldmMT^^Mrl^nipli using the IRAS PSCz survey were in good 
agreement with the predictions from the gravitational instability picture and provided new constraints on the bias for these 
samples. 

These studies have had much success at measuring the 2PCF and 3PCF of galaxies. Furthermore N-body simulations 
have been able to make precise predictions, given the cosmological parameters, of the 2PCF and 3PCF of the mass density. 
However until recently the galaxy-mass correlation functions have remained unmeasured. The measurements of the 3PCF 
of galaxies in fact can be used to constrain the large scale bias between galaxies and mass but the smaller scale bias has 
mostly remained a mystery. The fact that on small scales the observed 2PCF of galaxies and the 2PCF of mass from N-body 
simulations disagree leads to the requirement that there must be scale dependent bias between galaxies and mass if the CDM 
framework is correct. Another indicator of this complexity is that different populations of galaxies have different measured 
correlation functions and therefore are biased differently relative to each other. There is no reason why any particular sample 
of galaxies should trace the dark matter exactly. Trying to measure and model the bias using galaxy-mass correlations is one 
of the main goals of galaxy-galaxy lensing and the galaxy-mass 3PCF should help achieve this goal. 


3 MEASURING CORRELATION FUNCTIONS WITH WEAK LENSING 


The general relativistic lens equation relates the bending angle of light rays to the mass distribution. In the thin lens ap¬ 
proximation the deflection can be computed from the projected surface mass density Tj{x,y) = f dzp{x,y, z) where 2; de¬ 
notes the coordinate along the line of sight. It is customary make this dimensionless by dividing by the geometric factor 
^crit = cJDs/A'kGDlDls, where the distances are the angular diameter distances to the lens, to the source, and from the 
lens to the source. This defines a dimensionless 2D mass density, or convergence, K,{x,y) = T,(x,y)/'Ecrit- It is useful to define 
a 2D potential 

<j){x) = — j d^x'K,{x)\n\x — x\ (5) 

One can show that k and (j) satisfy Poisson’s equation 

2k = V ^ 4 > = 0,11 -I- 0,22 (6) 

where the subscripts with commas denote derivatives transverse to the line of sight. The other second derivatives of 0(5:) 
define the shear, 


271 = 0,11 - 0,22 and 72 = 0.12 (7) 

The two-component polar 7 is related to the anisotropic stretching of the source galaxy shapes in a linear way. With the 
ellipticity components defined in terms of the second moments of the galaxy light then it can be shown that {a) = 2"fi. Thus 
measurements of galaxy ellipticities provide an estimate of the shear. These measurements are noisy because galaxies are 
intrinsically elliptical with « 0.3. Theref ore accurate shear measurements require one to average over enough galaxies 

to reduce this “shape noise” to ((efj/A)^^^. See iBartelmann fc Schneider! J200lh for a review. 
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Because k, 71, and 72 are all related to (f) it is possible to determine k from the measurable 7. One such relation is 


Vk = 


7 1.1 + 72,1 

72.1 - 71.2 


( 8 ) 


To solve for k one needs to solve this partial differential equation, which states that the derivative of k is locally related to 
the derivative of the shear. By Fourier transforming this equation (or Equations | 7 | and one can see that k is locally related 
to the shear in Fourier space, but in real space it is related to the shear non-locally (i.e. it involves a convolution). There is 
also the problem of the mass sheet degeneracy: k can only be measured up to a constant with shear data. 

Wea k lensing first became a useful tool in cosmology as a method of measuring the projected mass around rich clusters of 

S ies llFahlman et al .ll994l : lciowe et al.ll99dljoffre et alJ200fll : llrgens et al.l2nndlTvson fc Fischeill99,‘Tl : lLuDDino fc Kais^ 

. Rich clusters are massive enough that the surface mass density is close to the critic al surface mass density 'Scrit 
(sometimes greater) and so there are sometimes giant arcs llHa.mmetlll99lt iKneib et al and even if not, the lensing 

shear is often very large and so can be measured reliably. Various methods of inverting th e shear map to produce 2D mass 
maps for clusters have been developed llLombardi fc BertirJll99^ : ISeitz fc Schneid^ l200lh . based on the original work by 
iKaiser fc Sauirej Jl99ijl . We will refer to any of these algorithms simply as KS algorithms. 


4 GALAXY-GALAXY LENSING 


Galaxy-galaxy lensing (GGL) denotes the lensing of background source galaxies by individual foreground lens galaxies. In 
this case, one does not have large enough signal-to-noise ratio to reliably measure the signal around each lens galaxy; in¬ 
stead one must a verage the signal f or many lens galaxies and be content with this statisti cal measure. The first attempt 
to measure GGL JTvsou et alj ll9j3) found only an u p per limit. The first d etection was bv iBrainerd et all (ll99(J) followed 
lTyj£)e11Ajij^io_^^;Tvsonr ir99^ : iGrifhths et alj ill 99fti ; iHiidson et all jl M^. The first very significant measurement was by 
EacheL£^n~il2QQl]l) in th e SDSS followed by o ther high S/N result s Iwil.son et alJlifKvIl : ISniith et aklbnoit llVlc.Ka.v et all 
^OCnTlHoekstoa^t alJl200^ . The latest result bv Isheldon et all (l2004l) makes high signal-to-noise ratio measurements of the 
3D galaxy-mass correlation function and the bias. 

The early GGL studies concentrated more on small scales {R < 1 h~^ Mpc) and were interested in constraining the 
mass to light ratios by modeling the mass profiles as isothermal spheres and extracting a velocity dispersion. However the 
best way to interpret GGL, especially at larger scales, is as a measurement of the galaxy-mass two-point correlation function 
(GM2PCF), ^gm{r)- The average density of mass around a galaxy can be written 

PmlG — T ^pm(r)] (9) 


After dropping the constant mass sheet which produces no lensing, the 2D mass over-density is just the projection 


E(R) — Pcrit — Pcrit Hn 


^iR) 


( 10 ) 


where 2 denotes the radial coordinate; i.e., F = B? + 7?. This defines the projected GM2PCF, Wgm{R), where R is the 
distance from galaxy to mass particle projected in the plane perpendicular to the line of sight. When the average mass profile 
is circularly symmetric, as in this case, the mean tangential shear 7 t in an annulus of radius R can be expressed in terms of 
E(i?) as 


^crit-fHR) = AE(R) = E(< R) - E(R) 


( 11 ) 


where E(< R) denotes the average E inside a circle of radius R and E(i^) is the mean in a narrow annulus of radius R. 
Written out in terms of Wgm, this is 

rR 


^E(/2) — Pcrit 




dR' R'Wgm{R') - Wgm{R) 


Differentiating with respect to R, one arrives at 


o 

Pcrit “ 


(R) 


dAE(R) ^ ^A^{R) 


( 12 ) 


(13) 


dR dR R 

One could then integrate this equation to get Wgm[R) in terms of tangential shear measurements. The unknown constant of 
integration is simply due to the usual mass sheet degeneracy. This equation is the circularly symmetric equivalent to Equation 
|H]and in fact can be derived by writing that equation in polar coordinates and assuming k is only a function of R. However 
rather than integrating this equation to get Wgm{R) one can calculate the 3D ^gm(r) directly . Since £„m.(r) is spherica lly 
symmetric, there exists an inversion formula for the projection called the Abel integral formula llBinnev fc Tremainelll98'^ . 


^gm{r) = — 

TT 


’ ^j^dWgrn 


(R) 


1 


dR y/R? — r'^ 


(14) 
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This equation was also used bv ISaunders et ^ dlQQSh to obtain the galaxy-galaxy auto-correlation function from the measured 
projected correlation function in the IRAS survey. Equation can be substituted into Equation LUJ to compute the gal axy 
mass correlation function from the shear measurements. This inversion method was recently used in isheldon et al.l (l2004l to 
make the first direct measurements of the 3D GM2PCF from weak lensing. 

One thing to note about Eniiation lMl is that one does not have data out to infinitely large R. If one has data from scales 
Rmin to Rmax then One can measure ^gm{r) only on these scales, and one can perform the integration only up to Rmax, 


^ Rmin) — 


-1 




(R) 


where C{r) is a correction 


C{r) = — 

7 T 




(R) 


dR _ J.2 


1 


+ C{r) 


dR 


JR? - r2 


(15) 


(16) 


To evaluate the correction term, one can either extrapolate Wgm{R) beyond the data region or use a model for ^gm{r). 
Either way, this results in a model dependence or uncertainty in the correction. Fortunately, because ^gm{r) drops off rapidly 
(typically like r~^'^ or faster at large scales) this correction term is negligible for an y reasonable £„m.(r) e xcept for r close 


to Rmax■ The first applicaion of th is inversion method t o galaxy-galaxy lensing was Isheldon et al.l (|2 
treatment of the method is given in l.lohn.ston et aJ ] (l2004ll . 


and a complete 


5 MEASURING THE HIGHER ORDER FUNCTIONS 

Just as the galaxy-mass two-point correlation function can be measured with the lensing signal around a galaxy, the galaxy- 
galaxy-mass three-point correlation function (GGM3PCF) can be measured with the lensing signal around pairs of galaxies. 
Here we derive the convergence field k around pairs of galaxies separated by a projected distance Ri 2 - 

In 3D, the probability of finding two galaxies and a mass particle in some triangular configuration in three infinitesimal 
volume elements, each of common size dV, is given by Equationl^ We call this dPggm{ri 2 ,r\ 3 ,r 23 ), where 1 and 2 indicate 
the galaxies and 3 indicates the mass particle. When we already have a pair of galaxies and just ask what is the probability 
of finding a mass particle, this is the conditional probability 

dP{m\gi,g 2 ) ^ dPggm/dPgg (17) 

Dividing this expression by dV and multiplying by the mass per particle yields the average density of mass around pairs of 
galaxies separated by ri2. 


pm (^13?'^231^12) — pcrit^n 


— pcrit^r) 


1 + ^39(^12) + ^gmjris) + ^gm(r 23 ) -|- Cg 

1 + 69 (’'12) 

^gm (fis) + C gm (^23) + C ggm 


l-b 


l+^ggiri 2 ) 


(18) 


With weak lensing, we find galaxy pairs at some projected distance R 12 and want to know the 2D projected mass density 
around these galaxies. So we need to project the mass over the radial direction but also average over the unknown radial 
distance between the two galaxies. Even though we may have measured redshifts for both galaxies, they cannot be used 
to cal culate this radial distance precisely because the radial peculiar velocities cause redshift distortions. Following [Pegb^ 
il98nll . where a similar projection is performed, we define u to be the unknown radial distance between the two galaxies and 
V to be the unknown radial distance between galaxy 1 and the mass particle (labeled 3). The radial distance between galaxy 
2 and the mass particle therefore is {v — u), since the three points define a triangle. With these two numbers u and v we can 
write the 3D distances in terms of the 2D projected distances 


2 ^2 , 2 

ri2 — Ri 2 + u 

f 13 = rtl 3 + w 

J'23 — R23 + (^^ “ w) 

To average over u, we need the normalized probability distribution P{u\Ri 2 ), 


P{u\Ri2 ) 


1 + Iasi's/ + ^ 12 ) 

J du{l+ £,gg{^JV? + RI 2 )) 

1 + C9g(\/^^ + -^^12) 

Lu -\- 'W\ 2 {Ri 2 ) 


(19) 


( 20 ) 


( 21 ) 
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One can see that this distribution is improper in the sense that one cannot normalize it over (— 00 , 00 ). Rather one needs to 
specify a maximum radial separation between galaxies to do the integration, which we will call Lu/2. We will similarly take 
the maximum range for the v integral as Ly 12 . 

Now we can write down the equation for the 2D mass density around galaxy pairs at separation Rn by taking Equation 
CHI and projecting this over v and averaging over u, 


S(Ri 3 , R23IR12) = pcrit^n 



pLy f2 


rLv/z 

— Pcrit^m 

1 

dv 

/ duP{ 

J 

-Ly/2 


^ — Lu 

— Pcrit^m 

Ly -I- 

Lu 

Lu , 

F- m 


+ Wgg(Rl2) 


1 + 


dv du P{u\Ri2) Pm(j' 13 ,r 23 |ri 2 
^gm (ns) +c gm (ns) + c ggm 


{Wgm{Rl 3 ) + Wgm{R 23 )) + 


'^ + ^gg{r 12) 
Z„ 


Lu + Wgg{Rl2) 


where we have defined the projected correlation functions, 
Wgm (Ri3 ) 


Wgm (-R 23 ) = 

Wgg{Rl 2 ) = 
Zggm{Rl 2 , R23, R23) = 


j dv + RI 3 ) 

j dv ^gmiV'^^ + ^^ 3 ) 

[ du ^gg{\/u^ RI 2 ) 

du dv (ggm{u,V, Ri2, R23, R23) 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


From here on we will ignore the constant term oc Ly in Equation CH since it is not observable with shear measurements. Since 
this is the only term which involves Ly explicitly, we can now let Ly —> 00 so the integrals involving v are now over all 
space. We can rewrite Equation 1221 without the constant as 

pcritLljn 


S(i?i 3 , R23IR12) = 


1 + 


W„g(Rl2) 


U)gm{Rl 3 ) + Wgm{R 23 ) + 


Lu - 


(27) 


For lensing we need the dimensionless mass density k, and so we divide by Ticrit- This Euclidean projection and division 
by Ticrit implictly assumes that the correlation functions do not extend over large enough distances to change the focusin g 
strength. A more general expression would include the radially varying PiFt inside the projection llXakada fc ,Tainll2003lJl . 
However, the correlation functions are only appreciable over several Mpcs, while typically varies over several hundred 

Mpcs and so this is a reasonable approximation. 

Given galaxy redshifts zi and Z2 and a source galaxy redshift Zs, Equation EH yields 


k(7?13, _R23|-^12) = 


Pcrit^n 


1 + 


Wgg(Rl2) 

Lu 




^crit (- 2^15 Zs') ^crit{z2^ Zs) Lu ^cr-it(.Z^ ^ Zs) 


(28) 


Here we have assumed that 21 ~ 22 in the third term, since Zggm vanishes if this is not the case. In most cases, we do not 
know the redshifts of the source galaxies precisely, so we need to integrate over their probability distribution P{zs) which can 
usually be estimated, e.g., using photometric redshift estimates. We define the univariate "Ecritizi) as 


GT — dZsPlZs)EcritlZi^ Zs) 

Now we can rewrite Equation 1281 as 


(29) 


K(i? 13 , i? 23 |^ 12 ) 


pcrit^n 


1 + 


Wgg{Rl 2 ) 

Lu 


'IVgmi^Rls') ^ 't^gm(-f^23) ^ Zggm 
^crit ('2l) EcruGT T[7 EcritGG 


(30) 


This is the main expression we will use to predict the shear around pairs of galaxies separated by i?i 2 . N-body simulations 
will be used to calculate the projected galaxy-mass correlation functions, Wgm, and Zggm, and Equation EOl then predicts 
k{Ri 3 , R 23 \Ri 2 ) around galaxy pairs for a survey of a given depth. 

Equations|Hlto[7|tell us how to take a convergence k and compute the resulting shear map. Furthermore we can apply the 
KS algorithm to the shear map to reconstruct the k map. Equation 1301 also tells us how to interpret this reconstructed k map 
in terms of projected correlation functions. Since the term Wgg{Ri 2 ) is just the projection of the galaxy-galaxy correlation 
function ^gg, which we can measure, and since we can also will measure Wgm by the usual two-point galaxy-galaxy lensing 
method (described in Section |1J , we can measure the 2D three-point function Zggm directly from these other measured 
functions and the reconstructed k. 
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If one lets the maximum radial galaxy separation, Lu, go to infinity, one does not get any contribution from Zggm- This 
just says that with an infinite background of galaxies, the ratio of physical galaxy pairs to random projected galaxy pairs 
will go to zero. If one’s goal is to measure the GGM3PCF, then one should not make this measurement around pairs with 
widely different redshifts: they only give contributions from their projected two-point terms Wgm{Rii) and Wgm{R 2 i), which 
provides no new information beyond what one has already measured with two-point galaxy-galaxy leasing. 

Instead one should set Lu no larger than the scale of measurable galaxy correlations. A good value might be Mpc, 

since the 2PCFs are small beyond that and dominate over the 3PCF though one might want to make Lu larger when making 
measurements at larger scales. At a minimum, one will want Lu greater than a few times the maximum projected scales at 
which one is making the measurements. In fact one could choose Lu smaller to attempt to maximize the signal-to-noise ratio 
of the measurement. However at smaller scales redshift distortions make the interpretation of radial distance problematic and 
so this probably will not prove useful. Making cuts on redshift differences at any scale should really be interpreted as a “fuzzy 
selection” in physical space. This selection function <?!>(u) can be modeled with N-body simulations and incorporated into the 
integrals above, rather than assume the top-hat form, to better interpret the measurements. 

We can now estimate at which scales the three-point function Zggm will dominate over the two-point terms in Equation 
EH For simplicity we choose the isosceles triangle configuration so that R \2 = Ria = R 23 = R- Let us model the three 2PGFs 
as power laws, ^gg{r) = ^gm{r) = (r/ro) *, ignoring bias parameters which are of order unity. We will take ro = Qh~^ Mpc 
and 7 = 1.8. The three projected 2PGFs are all then given by w{R) = 3.7ro(ro/R)~''~^. We will further assume that the 
projected three-point function Zggm can be written as Zggm = Qlwgg(R)wgm(R) + Wgg(R)wgm(R) + WgmiRJwgmiR)], with Q 
approximately constant. From Equation 1271 we can write the three-point mass density S(R) for this configuration as 


S(R) = 


Pcrit^n 


I 3^7_ro /lO ^^7-1 
Lu ^ -R ^ 


3.7 ro 


(s) 


7-1 


2-b3g 


3.7 ro /ro\T'-i 


(S)' 


(31) 


Since the usual two-point contribution from each galaxy, which we will write here as Tj 2 {R), is Tj 2 {R) = PcritLlm 3.7 ro(ro/i?)^ ^ 
lEauation lKH . we can write Equation 03 as 


E(i?) = 


S2(i?) 


1 I 3^7_ro ,'10 W-l 
^ ^ Lu ^ R ^ L 


2 + 3Q 


3.7 ro /roX-^-i 


m 


(32) 


The first term in brackets comes from the two-point terms in Equation 03 The next term is due to the three-point function 
and equals the first term at R = Rug, where 


Reg — ro 


2L, 


3.7 3 Qro 


0.9h"^Mpc 


(33) 


where we have taken Q = 1 and Lu = 30/i“^Mpc. For small scales R <C Reg, where the three-point term dominates, we see 
that E(i?) ~ 3 Q ^ 2 ( 7 ?), independent of Lu- In general we want to subtract off the two-point term so we define the reduced 
three-point surface mass density 

E., 

We see that at R = lQh~^ Mpc, T,^{R) = 0.5 T, 2 {R) and at R = 100h~^ Mpc, ^ 3 ( 7 ?) = 0.1 E2(7?). Since we can measure the 
GM2PGF to at least R = 10h~^ Mpc in the SDSS Jsheldon et al.ll2004ll . we should be able to measure the GGM3PGF if we 
have a similar number of galaxy-pairs as we have galaxies. We will study this further in Section EOl 


6 HALO OCCUPATION MODELS 

The halo occupation distribution models (HOD) p rovide a simplified prescription for populating simulated or analytic dark 
matter halos with galaxies. In the HOD framework llBerlind fc Weinbergl20o3lBerlind et al.l2003il . one makes the assumption 
that the number of galaxies, N, that occupy a dark matter halo is drawn from a distribution P{N\M) that only depends 
on M, the mass of the halo. A priori, one might imagine that this is too restrictive and that the probability distribution 
could depend on other parameters such as the shape of the halo, the formation time, the local density of other halos, or a 
host of other parameters. However, hydrodynamical simulations indicate that halo mass is the controlling factor. Once one 
has determined the number of galaxies in each halo, one needs to decide on how to distribute the galaxies within each halo. 
Together these two choices define a complete halo occupation prescription. 

The chief usefulness of the HOD model for large scale structure studies is that it provides a model for galaxy bias and 
therefore allows one to compute the various galaxy correlation functio ns in a straightforward m anner. This technique was first 
used with ra ndomly placed clus t ers w it h specifi e d radia l profiles bv iNevman fc ScottI ^1952^ to predict clustering statistics 
and also by iMcClelland fc silS (Il977ll : IPeeblesI Jl974lJl used similar models to compute the galaxy correlation function. 
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More modern treatments benefi t from knowledge of the halo mass fu nction Jpress fc Schecht ^ll974l: Isheth fc Tormerilll99fll : 
I.Tenkins et al.ll200ll : lwhi't3l2002h and halo mass profiles JNavarro. Frenk. fc Whit^^Q^ that are seen in N-bo dy simulations. 
Another reason for the recent revival of this technique is the success of semi- analytic models of galaxy formation ilWhite fc ReesI 


). These models, which try to simulate the complex 
physics of galaxy formation, produce HODs which when combined with N-bodv dark matter simulations, resu lt in a galaxy 


iKanffmanu et al 




Due to this 


auto-correlation function that resembles the observed power law 
improved knowledge of the mass function, halos profiles, and the HOD, th ere has been much recent wor k on computing the 
various N-point galaxy correlation functions with the halo model forma lism (ISeliakl200d : lMa fc Frvl2nnnal l3 IPeacock fc SmithI 
I 2 OO 0 I: Iscoccimarro et alJl200l'5 iTakada fc .TairJf2003al : IScrantonll200ijl . 

In the purely analytic halo model formalism, one combines the HOD with an analytic halo mass function and analytic 
halos profiles that are fit from N-body simulations. Another approach is to take an N-body dark matter simulation, identify 
halos, and populate these halos directly with galaxies using the HOD prescription. This is the method we will use in this 
study. This method yields a set of dark matter positions and a set of galaxy positions that one can use to calculate the various 
correlation and cross correlation functions. It has the advantage that it samples a realistic distribution of halo shapes and 
structures not captured by the analytic method, which usually assumes spherical halos with a universal profile. In particular, 
while the constraint of spherical halos should not matter for the 2PCFs, it may be more important for the higher N-point 
functions since they are anisotropic. In fact, the analytic halo models do not agree precisely with N-body predictions for the 
3PCF, because the former are not hierarchical on small scales. 

However, populating dark matter simulations with galaxies also has its drawbacks. One problem is that if one wants to 
change the dark matter power spectrum or vary other cosmological parameters, one needs to run another N-body simulation. 
With the purely analytic technique, one can calculate the various correlation functions simply with numerical integrations; to 
change cosmological parameters one just has to re-run these integration routines, which take far less CPU time. Thus with the 
purely analytic technique one can explore parameter space more efficiently. Also with the purely analytic technique one can 
compute the various statistics with no errors. With the N-body population technique one always has shot noise and sample 
variance. However if one has a large enough volume and enough particles and galaxies, one can make predictions that will have 
theoretical errors smaller than errors coming from real survey data. Another thing to note about the analytic models is that 
the 2PCFs are easily computed with simple one dimensional numerical integrations but the higher order N-point functions can 
require much higher order numerical integration. For example, to compute the full three-point function one must compute a 
seven-dimensional num erical integral and since this is generally intractable, one has to resort to some approximation methods 
llTakada fc ,IairJl2003a^ . In addition to the extra numerical computation, the mathematical expressions become significantly 
more complicated and require more assumptions about parameters and input from simulations. 


7 THE N-BODY SIMULATION 


The N-body dark matter simulation used for this study is the public simulation made available at http://pad.berkeley.edu/Siml/ 
This simulation has 51 2^ particles in a pe riodic cube of length 300 h~^ Mpc. The simulation is evolved from z = 60 to z = 0 
using a TreePM code llWhitell20n 3. l2003h . The initial conditions are set using the Zel’dovich approximation to displace the 
particles from the uniform grid. The cosmological par ameters used are Qm = 0.3, Q,a ~ 0.7, h = 0.7, n = 1, Q.b h? — 0.02 
and (Tg = 1. The CDM transfer function used is that of lEi.senstein fc m] (11991^1 with the baryon wiggles smoothed over. The 
simulation has an effective Plummer force-softening scale of 20 h~^ kpc which is fixed in comoving coordinates. The mass of 
each dark matter particle is 1.7 x 10 ^°/i“^Mq. The simulation outputs three positions and three velocities for each particle 
but only the posit ions are used for this study. The dark matter halos are found using a Friends-of-Friends (FoF) algorithm 
IDavis et al.lll98^ with a linking length of 0.2 in units of the mean inter-particle separation. 

The N-bodv simulation is then populated with galaxies at « = 0 using halo occupation distribution (HOD) techniques 
based on iRerlind fc WeiubergI (l2002ti . courtesy of Berli nd. The HOD param eters are chosen to produce a galaxy 2PCF that 
matches ^gg{r) of the SDSS volume-limited sample of lZeha.vi et alJ (1200,31! for absolute magnitudes Mr < —21. There are 
123081 galaxies corresponding to a comoving number density of Ug = 4.56 x lQ~^h ^/ Mpc? . This number density corresponds 
to an absolute magnitude cut Mr < —20.2 according to the luminosity function of iBlanton et al.l (|20oJ. The P{N\M) used 
is one that has a mean {N{M)) which is as follows 


(7V(M)) = 



M < 

M^ir^ < M < Ml 
M > Ml 


(34) 


The parameters we use are Mmin = 1-09 x 10^^h~^M q and Mi = 9.97 x 10^^h“^MQ, with /3 = 0.83. 

Equation 1341 gives the mean number of galaxies for a halo of mass M. We must allow some scatter when we choose the 
number of galaxies for each halo. Since we must choose an integer number, we choose either one of the integers bracketing the 
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mean with relative probabilities that keep the mean constrained to equal Equation 01 when the mean is non-integer. When 
(N) is an integer, we simply choose that number with no variance. This leads to a variance that is substantially sub-Poisson 
and in fact is the minimum variance that is not identically zero. To be more specihc, if the mean galaxy number is N = N + e 
with N an integer and e the fractional part, then we choose N with probability 1 — e and -|- 1 with probability e. One can 
easily show that the resulting variance is = e(l — e). The variance therefore depends on where between two integers N 
happens to lie, and it is maximized at — 1/4, half way between two integers (e = 0.5). 

Once the number of galaxies is chosen for a given halo, the hrst galaxy is always placed at the halo center of mass and 
the remaining galaxies are each pla ced randomly at the positions of one of the othe r dark matter particles in that halo. For 
more details of this HOD model see lZehavi et al.| (l200,‘^l : lBerlind fc Weinberel ll2nn2h . 

8 CORRELATION FUNCTIONS FROM THE SIMULATION 
8.1 The 3D two-point functions 

We measure the 3D two-point correlation functions in the simulation using two different methods. On large scales we bin 
the data onto a cube of 256® pixels, Fourier transform to measure the power spectrum, and then transform back to get the 
correlation function. This method works well but because of the finite resolution imposed by memory limitations, it is limited 
to measuring only on large scales. This method makes reliable measurements of ^(r) on scales 2h~^ Mpc to 50/i“^ Mpc. 

For smaller scales, we do simple pair counting to estimate ^{r). We break the entire region of data into A®gg different 
segments and write each to a different file. This enables one to search quickly for pairs since one only needs to read in the 
neighboring 27 (or 9 in 2D) segments to measure pairs within a separation of 1 segment. We typically use Weg = 100, which 
gives a maximum search radius of 300/100=3/i“^ Mpc. This range overlaps the range of the Fourier method, and we hnd 
excellent agreement between the two methods in the overlap region. 

The 3D correlation functions are shown in Figure^ These plots show that on scales greater than about lh~^ Mpc, the 
three curves trace each other fairly well. However, on smaller scales both ^gmir) and ^gg(r) are biased positively with respect 
to the mass. The bias of ^gg(r) is expected, since it is neccesary to have substantial bias at small scales to fit the observed 
^gg(r), which is known to be close to a power law. Also one can see that on these small scales, ^gm(r) is greater than ^gg{r). 
Both of these qualitative results are expected with the HO D framework. On small scales, the 2PCFs are dominated by the 
contribution of galaxies in the smaller halos llSeliakll200(ill . which are the most abundant; many of them contain only one 
galaxy at the center of the halo. Because the galaxy is at the center, ^gm{r) traces the cusp of the dark matter halo and will 
therefore be steepest at small scales. ^mm{r) is smoothed out somewhat, because this is basically the halo profile convolved 
with itself. In the case of ^gg(r), it is suppressed at small scales due to the fact that small halos have a P(N\M) that excludes 
a second galaxy, so ^gg only gets contributions from bigger halos and the two-halo term (galaxies in different halos). This 
argument explains qualitatively why ^gm{r) is the largest at small scales. The most likely reason why ^gg{r) is larger than 
is that the P(A|M) gives more weight to smaller halos for ^gg(r) and one galaxy is always at the center. These effects 
must be more important than the exclusion effect described above, for ^gg{r) to be larger than 


8.2 Scale dependent bias 

It is conventional to define the bias parameters as ratios of the correlation functions or power spectra. The bias parameters 
can therefore be scale dependent; in fact scale dependent bias is necessary if one is to get a nearly power-law galaxy corre¬ 
lation function (that one sees in galaxy surveys), from a non-power-law mass correlation function that one sees in N-body 
simulations. We dehne the gala xy bias pa rameter as b^{r) = ^gg{r )and the galaxy-mass cross correlation coefficient 
r(r) = igm.{r) j JPenlllQQ^ . From these definitions one can also see that 6(r)/r(r) = iggir)/^gm{r). The ratio 
b/r is directly measu rable from combining galaxy-galaxy leasing measurements with galaxy auto-correlation measurements 
llSheldon et aljUoOdll . 

These three quantities 6, r, b/x for the simulation are shown in Figure]^ The bias parameter b starts at about 1.5 at 
r = 0.05h“^ Mpc, decreases to about 0.95 at about r = Mpc and then increases slightly at larger separation. This 

behavior is most likely due to the transition from the one-halo term to the two-halo term. At large scales, b plateaus at about 
1.05, though there is some indication that it falls off slowly toward larger scales. The r parameter is decreasing monotonically 
with r out to 2 h~^ Mpc and beyond that it is flat at r = 1. 

The ratio b/r is shown as the dot-dashed curve. It is fairly flat from r — 0.2 h~^ Mpc to 1.0 h~^ Mpc and then rises to a 
slightly higher plateau by 3.0h~^Mpc, beyond r = 2 h~^ Mpc b/x — b since r = 1. 

On small scales r < 0.2h~^ Mpc, b/x drops by almost half. This scale dependence of b/x would not yet be apparent 
in SDSS lensing data but it should become measurable i n future surveys. The overall amplitude b/x « 1 agrees with SDSS 
galaxy-galaxy lensing measurements llSheldon et ahlEoPlIl for Qrn ~ 0.3. 
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Figure 1. The three 3D two-point correlations functions from 0.06/i ^Mpc to 50 h ^Mpc in the N-body simulation. The solid curve 
denotes fgm(T), the triangles denote ^gg{r), and the squares denote 




Figure 2. The bias parameters as a function of scale. The solid curve is b = y^^ggl^mm- The dashed curve is r = ^gml^gg^mm and 
the dot-dashed curve is 6/r. 


8.3 The 2D two-point functions 

We also measure the 2D two-point functions where we first project the data (i.e. ignore the third coodinate) and measure 
correlation functions of the 2D densities. However these are simply projected versions of the 3D two-point functions and the 
analytic projections of the 3D functions agree very well with the measured 2D functions. We note that projected correlation 
functions, which are usually referred to as Wp, have dimensions of length. If one measures a 2D correlation function by first 
projecting the points to 2D and then computing the 2D correlation function W 2 d in the usual way, the result is dimensionless 
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Figure 3. The three projected two-point correlations functions from R = 0.05 h '^Mpc to 10 h '^Mpc. The solid curve denotes Wgm{R), 
triangles denote Wgg[R), and the squares denote Wmm{R)- 


and it is equal to Wp/L, where L is the dimension of the box in the radial direction. When we refer simply to w, we mean Wp, 
which has dimensions of length. Figurej^shows the three projected correlation functions. 


8.4 The 3D three-point functions 


Next we measure the 3D 3PCFs (ggg and Cggm- Our main purpose for measuring these functions is simply to check our 
results against other published results of Cggg- For the most part we will be interested in the 2D three-point functions, since 
they are directly measurable with leasing. The three-point function depends on the three triangle sides, however it is usually 
parametrized with the three variables being two side lengths r and q and the angle ip the intersection of these two sides. 
Motivated by the hierarchical form of Equation |11 it is customary to divide Cggm by the usual symmetric combinations of 
2PCFs to define the scaled 3PCF Qggm for two galaxies and a mass particle 

Cggrri 


^ggjl2)^gm(ri3) + ^gg{l'12)^gm{'r23) + £,gm{l'13)P,gm{l'23) 
and likewise for three galaxies 

Cggg 


Qggg — 


^99 (ri2)C 99 (ns) + ^99 (ri2)C 

99 {r23)+^ 99 {ri3)i 99 J 23 ) 


(35) 


(36) 


We compute the 3D 3PCFs only on small scales by first creating a galaxy pair catalog; for each pair, we search for 
neighboring galaxies or dark matter particles to complete the triangle. We limit the first galaxy pair separation r from 0.1 
h~^ Mpc to 0.9 h~^ Mpc and search around the center of the pair to a maximum radius of 3.0 h~^ Mpc. 

Figure 2] shows the reduced function Qggm for galaxies with separations r = 0.25 ± 0.04/i“^ Mpc for 5 side ratios q/r 
from 1 to 5. Figure|Hlis the same plot but for galaxies with separations r = 0.5 ±0.04h“^ Mpc, and FigureElfor galaxies with 
separations r — 0.9 ± 0.04h“^ Mpc. 

We also plot the function Qggg, although this is much noisier since there are many fewer galaxies than dark matter 
particles. These three plots are Figures 13 H and|3for galaxy separations 0.25, 0.5, and 0.9 h ^ Mpc respectively. 

Although there is a lot of variation in these plots, there is a generic trend that describes all of them. First, the hierarchical 
scaling works fairly well, in that both Qggm and Qggg only vary between 1 and 3 over all of the plotted ranges. Except for 
the case q/r = 1, the amplitude of Q is higher on either end when ip = Q ox ip = n and lower in the middle when ip = 7r/2 
and is fairly symmetric about ip = -k 12. This indicates a small preference for collinear triangles. This preference for collinear 
structures is pred icted in perturbation theory and observed in galaxy s urveys for much larger scales but with a much larger 
variation with ip iFrvlll98'3 : iTakada fc ,Tainlf2003at iFeldman et al.ll200lh . Halo model calculations by iTakada fc .TainI ll2003^ 
(hereafter TJ) reveal a flatter curve with ip for their Qggg then ours at r = 0.1h“^ Mpc. This small difference may be an 
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Figure 4. The 3D Qggm as a function of angle -i/; for galaxy separation ri 2 = 0.25h ^ Mpc and for various q/r ratios. 
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Figure 5. The 3D Qggm as a function of angle i/i with galaxy separations ri 2 = 0.5h ^ Mpc and for various q/r ratios. 

indication of halo ellipticity, since TJ are constrained to use spherical halos in their analytic calculation, or it may simply 
be due to TJ using different parameters from the ones in this simulation. The larger scale r = lh~^ Mpc plots of TJ shows 
more complicated behavior where the curves for different q/r are not self similar. This is due to the fact that all three terms , 
one-halo, two-halo, and three-halo, are contributing at this scale and each term contributes differently for each q/r. Our curves 
at r = 0.9h,”^ Mpc look different from theirs and not too different from our own at smaller scales. This may indicate that 
in our simulations the one-halo term still dominates at this scale whereas in TJ the two and three halo terms are becoming 
important. 

We also see that the q/r > 1 curves are fairly self-similar and simply decrease in amplitude as q/r increases. This is also 
seen in TJ at small scales. For the case q/r = 1, we have somewhat different behavior: the curve is not symmetric about 
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Figure 6. 


The 3D Qggm as a function of angle ip with galaxy separations ri 2 = 0.9h ^ Mpc and for various q/r ratios. 
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Figure 7. The 3D Qggg as a function of angle ip with galaxy separations ri 2 = 0.25h ^ Mpc and for various q/r ratios. 


Ip = -E12. For large angles {ip ~ tt) the curves looks similar to the other q/r > 1 curves but at ip close to zero Q decreases 
substantially. These t/) ~ 0 triangles are such that the third galaxy or dark matter particle is right on top of the first galaxy. 
So this decrease in Q indicates that the hierarchical scaling does not work as well very close to the two galaxies. The value 
of Q can get as small as 0.1 for cases where the third particle is very close to the first (or second) galaxy, but it never goes 
negativ e at the scales we consider. It is possible in principle for Q to be negative; for example, in TJ and also in perturbation 
theory i Frieman fc Oa.ztanagairiQOljl . Q can go as low as -4 on large scales ( > 10h~^ Mpc). This decrease in Q at i/) ~ 0 is 
also seen in TJ and seems to be characteristic of the one-halo term and probably is affected by the inner halo slope and the 
number of small mass halos. 
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Figure 8. The 3D Qggg as a function of angle ijj with galaxy separations ri 2 = O.bh ^ Mpc and for various q/r ratios. 
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Figure 9. The 3D Qggg as a function of angle i/i with galaxy separations r = 0.9h ^ Mpc and for various q/r ratios. 


8.5 The 2D three-point functions 

The 2D GGM3PCF is the object that can be measured directly with weak lensing. As we described in Section|^the 2D mass 
over density, T,{x,y), around pairs of galaxies separated by a fixed projected distance i?i2 is given by Equation I3III which we 
rewrite as 


T.{x,y) = 


peril 


! + ■ 


g(-^12) 


Lu 




Zggm{x, y) 


(37) 
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where 

Ri 3 {x,y) 

= -\/(x + Ri2/2)2 -b t/2 

(38) 

R 23 {x,y) 

= \/(a; - I?12/2)2 -b t/2 

(39) 


Here we have simply used the parameter R12 for the galaxy separation and x, y for the Cartesian coordinates centered on the 
midpoint of the galaxy pair and rotated to make the galaxy pair lie along the x-axis. These three parameters R\2, x, and y 
specify a triangle connecting the two galaxies and a mass particle. 

We display contour plots for the 2D density around pairs with R \2 = 0.9h~^ Mpc in Figure irni These boxes are 4.24 
h~^ Mpc on a side. This is a projection through the whole volume for all projected galaxy pairs, so here Lu ~ 

Mpc. The top left contour plot shows T.{x,y). The top right plot shows the contribution coming from the two-point functions 
Wgm{Ri 3 ix,y)) -b Wgm{R 23 {x,y)). This must be subtracted from E(a:,y) to get the reduced three-point function. One can 
see how these terms are each radially symmetric about the two galaxies, so their sum creates most of the bimodal pattern 
seen in the first plot. The plot in the lower left is Zggm{x,y). Even after subtracting off the W13 + W23 parts, there is still a 
bimodal pattern which is expected from the hierarchical model. Finally the plot in the lower right shows the scaled three-point 
function, qggm{x,y), which is defined 


qggm{x, y) 


_ ^ggrnjx, y) _ 

Wgg(^Rl2^ 'Wgni(^Rl3^ ~\~Wgg(^Rl2^ Wgm(R23^ '^gm(^Rl3^ '^g7n(^R23^ 


(40) 


Here, the distances R\3 and R23 are functions of x and y fEquations 1331 and but Ri2 is fixed. At larger scales qggm, is 

fairly constant at 1.5 although it drops off slowly. On small scales around each galaxy, there are holes where qgg-m drops from 
~ 1.8 on the outside to ~ 0.6 at the center. This is similar to what we saw in 3D: the hierarchical scaling does not appear to 
work well very close to each galaxy. This may be a feature of the HOD models, since they do not associate sub-halos to each 
galaxy: non-central galaxies do not have a mass excess around them at small scales as do the central galaxies. In reality, it is 
probably the case that non-central galaxies have small sub-halos of dark matter around them. We do know at the very least 
that they have baryonic mass tightly clustered around them in the form of stars, so the HOD models without sub-halos must 
give incorrect galaxy-mass correlation function at these small scales (a few kpc) where baryons dominate. 

Figure rTTI shows the amplitude of qggm along the y ~ 0 axis, now for galaxies separated by R 12 = 0.25 h~^ Mpc. From 
this plot one can see that on larger scales qggm drops off slowly, but on small scales there are holes around the galaxy pairs. 
On scales larger than the galaxy separation, qggm is mostly radially symmetric about the center. In Figure m we show the 
radial profile of qggm about the center of two galaxies with R12 = 0.05/i“^ Mpc. On scales comparable to R12, Qggm is rising 
as it climbs out of the hole and at larger scales it drops off. 


9 EXTRACTING INFORMATION FROM SHEAR MEASUREMENTS 

As we saw in Section |11 we can use the method of galaxy-galaxy lensing (GGL) to measure the galaxy-mass correlation 
function. By measuring the tangential shear around lens galaxies, one can reconstruct Wgm{R) and then invert the projection 
to obtain ^gm{r). In Section |3 we calculated the 2D mass density around projected pairs of galaxies and showed how it was 
related to the GGM3PGF. We will now discuss how to measure this statistic from shear data using the SDSS as our fiducial 
data set. We assume we have a spectroscopic sample of galaxies with measured redshifts that we will use as our lens galaxies 
and a fainter photometric sample of galaxies without redshifts (though possibly with photometric redshifts) that we will use 
as source galaxies. 

Our first task is to find pairs of lens galaxies out to some maximum projected radius. We will also want to impose a 
cut on their redshift difference. As we discussed in Section]^ this cut imposes an effective cut in the radial direction Lu. 
One will probably want to keep this fairly large, such as 30 — 50 h~^ Mpc, to limit the effects of redshift distortions on the 
interpretation of the results. 

Once we have a galaxy pair, we rotate the local coordinate frame to align them along the x-axis. There was no need to 
do this with GGL since the average tangential shear around single galaxies must be isotropic. In that case, we just chose a 
set of radial bins and averaged the tangential shear for objects in each bin. Now we have a pair of galaxies which breaks this 
spherical symmetry and therefore need to deal with both components of the shear and keep track of both x, y components of 
the source galaxies. Thus, we now have a 2D problem and will measure a 2D shear map rather than just a ID shear profile. We 
choose an x,y grid on which to store the ellipticity measurements and use the rotated x,y coordinates of the source galaxies 
to calculate which pixel each source galaxy belongs to (for a given lens galaxy pair). 

We want to build up a map of the average source galaxy ellipticity stacked (averaged) over many lens galaxy pairs. 
However, since the lens galaxy pairs have a range of projected separations, we also need to bin in this separation, R12, as 
well and therefore produce a shear map for each R 12 bin. If we bin too coarsely, this will have the effect of blurring together 
different R 12 correlation functions and we lose information. If we bin too finely, the number of lens pairs per bin becomes 
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Figure 10. Top left: the 2D projected mass density Y,{x,y) for fixed projected galaxy separation R 12 = 0.9h~^ Mpc. At small scales 
there is a bimodal distribution with peaks around the locations of the two galaxies. At larger scales the density has elliptical contours. 
Top right: the contribution from the two point functions 1^13 -\-W 23 - These are simply computed from tabulated values for these functions 
and so this has no apparent noise. Bottom left: the galaxy-galaxy-mass three point function Zggm- This is obtained from taking Yi(x^y) 
and subtracting off the constant and the two-point contributions (see Equation 1 ^ . Lower right: the reduced projected three point 
function qggm, which is obtained by dividing Zggm by the symmetrized combination of products of two-point functions (Equation |42J. 
At large scales Qggm drops off slowly and retains some ellipticity. At smaller scales it reveals holes around the galaxies where qggm drops 
off sharply. These features of the halo models might not be present in reality since these models do not include sub halos of dark matter 
or even baryonic mass around the galaxy centers. All of these densities are log-scaled except qggm and the length of these boxes are 
4.24 h~^ Mpc on a side. 


smaller and the shot noise will overwhelm the signal. With real data, one would also want to bin the len galaxy sample by 
luminosity and perhaps by color, but we no not have that information in our simulation. 

The fact that this is a 2D problem with no exact symmetries means that instead of just a few radial bins one might have 
hundreds of pixels. In addition, as noted above, one has to separate the lens galaxy pairs into several bins of Ri 2 - Because of 
this, one cannot expect a high signal-to-noise ratio in any particular shear map pixel since the signal-to-noise ratio is only of 
order 1-10 in present measurements of the GGL signal from the GM2PGF with only a few radial bins. However, as we will see, 
one can use approximate symmetries to make measurements with signal-to-noise similar to the simpler GGL measurements 
of the GM2PCF. 

To calculate the shear map from T.{x,y) we solve for the 2D lensing potential 4>(x,y) by solving Poisson’s equation (see 
Section ??). Then we numerically differentiate <f) to obtain the shear components 71 and 72 . We will use as our model the 
density Tj{x,y) for R 12 = 0.8 h~^ Mpc. For simplicity, in calculating Ecrit, we choose all the lens redshifts to be zl = 0.1 and 
all the source redshifts to be zs = 0.3, similar to the mean lens and source redshifts used in the SDSS lensing measurements. 
We display the shear map for a box of width 6 h~^ Mpc calculated from this density in Figure ESI In Figure d we show a 
higher resolution image of the region close to the galaxy pair. This smaller region is 1.2 h~^ Mpc on a side. From the first 
plot one can see that at larger radius the shear is mostly tangential though one can notice a residual ellipticity to the pattern 
that drops off as one goes out further. 

In the central region, the shear map is more complicated and it helps to look at the second plot. Figure ITU To understand 
this pattern one has to remember that the shear is a vector (technically a spin -2 held) and this pattern is the superposition of 
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R (Mpc/h) 


Figure 11. A slice through the 2D qggm{x,y) along the y = 0 axis which is the direction passing through the two galaxy centers. The 
galaxy separation is R 12 = 0.25 h~^ Mpc. At small scales one can clearly see the holes around the galaxies (which may be an unrealistic 
feature of halo models). On larger scales qggm drops off slowly. 
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Figure 12. The radial profile of the 2D qggm for closer galaxy pairs, R 12 = 50 h~^ kpc. The inner part of the curve is still rising 
as qggm climbs out of the holes. At scales of a few times R 12 the profile is fairly radially symmetric and decreases slowly. There is an 
apparent bump at 1.5 h~^ Mpc where the two-halo term typically begins to dominate. 


two terms which are mostly tangential about the two galaxies. Because it is a vector, there can be constructive and destructive 
interference. The region in the middle is a region of constructive interference where the shear from each galaxy has the same 
phase (i.e., direction), so this is where the shear is the largest. There is also a ring of destructive interference where the shear 
from each galaxy cancels out the other and no net shear is produced. 

In Figure IT^ we show the individual shear components around the same 2D density map. Here the scale is 4.2 h~^ Mpc 
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Figure 13. The mean shear map around pairs of galaxies of separation R \2 = 0.8 h~^ Mpc. The size of the region is 6 h~^ Mpc on a 
side. At larger scales the signal is mostly tangential but in the middle the signal is more complicated. See next figure for a more detailed 
view. 


on a side. The top left shows k = 'SfT,crit- Top right shows the magnitude of the shear I7I = 7/71 +72- The bottom panels 
show 71 and 72 from left to right respectively. 

The K plot shows the bimodal pattern on small scales and elliptical contours on larger scales. The plot of I7I shows a 
dark ring where the net shear is zero due to the interference and the bright spots where the shear is largest in the center and 
just on the outsides of where the galaxies are. It also shows how I7I has elliptical contours on larger scales. 

In Figure ITCl we show the tangential shear 7 t and the radial shear 'yu defined as yr = —71 cos 20 — 72 sin 20 and 
7 _r = —72 cos 20 + 71 sin 20 where 0 is the angle from the midpoint of the x-axis. The top left is yx and top right yu. The 
bottom left is |7 t|/| 7| and bottom right is |7 t|/| 7|. In the bottom color scale, black is zero and white is 1. The bottom plots 
show how at large scales the shear is almost entirely tangential whereas at small scales the pattern is complicated and both 
7 t and yu make contributions. 

This is an example of how there is an approximate symmetry that allows one to ignore the radial shear component 
at large scales. On these scales one can just measure the tangential component and see how its magnitude changes as a 
function of R and 0. For example Figiire IT tI shows the tangential shear in an annulus (centered between the galaxies) of radius 
R = 2.0±0.1 h~^ Mpc, as a function of 0 (measured from the x-axis), divided by the mean value in the annulus. A quadrupole 
pattern is evident and we over-plot the function 7 t/7to( 0) = (1 -b e cos 20) with e = 0.07. One can relate the moments of t he 
shear map to moments of k through methods developed bv ISchneider fc BartelmanrJ (IiQQtII : iNataraian fc R.efregieJ J200fil . 


10 ESTIMATE OF SIGNAL-TO-NOISE RATIOS 

So far we have calculated the shear signal we expect around foreground galaxy pairs, and now we investigate the amount of 
noise we expect. The noise in a weak leasing measurement is almost entirely due to shot noise from galaxy sh apes, although 
on the largest scales, sample variance can increase the variance and add some covariance llSheldon et alJl20o3l . We will only 
discuss the shot noise since it is the dominant source of noise on scales less than a few Mpc. The shot noise on a particular 
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Figure 14. A zoom of the previous figure of the central region near the two galaxies. The scale is 1.2 h~^ Mpc on a side. In the middle, 
the mostly-tangential shear from each galaxy adds constructively. There is also a ring of destructive interference where the shear from 
each galaxy is in the opposite direction of the other. 


measurement of -ji is about 0.2/^/N, where 0.2 represents the RMS intrinsic ellipticity of galaxies divided by 2. The number 
N for GGL is the number of lens-source galaxy pairs in a given radial bin. For our analysis, we have pairs of lens galaxies and 
so N is the number of lens-lens-source triplets in a given measurement bin. We now estimate what this number N will be. 

If we have Np foreground galaxy pairs in a given projected separation bin, R 12 , and a 2D number density as of source 
galaxies, then the number of triplets is just N = NpOsdA where dA is some surface area element in which we estimate the 
shear. We have assumed that the source density as is uncorrelated with the density of foreground pairs, which should be 
correct if lens and source galaxies are at different redshifts. To estimate the number of lens galaxy pairs, we calculate the 
average number of neighbors per galaxy, (rip), in a thin annulus a bout projected separation R of width /S.R. This is just an 
integral over the galaxy-galaxy correlation function JPeebleslll98fll . 

_ 

(rip) = Ug 2^ R AR / du (1 + ^gg{\/ v? + R^)) (41) 

^—^xj 1^ 

= Ug 2 'k R AR {Lu +'Wgg{R)) (42) 


where Ug is the 3D number density of lens galaxies. The total number of pairs in this bin is then Np = NL{np) /2, where Nl 
is the total number of lens galaxies. The number Nl is just the number density rig times the comoving volume probed by the 
survey, V. Finally we can write the number of lens-lens-source triplets, N, in surface area element dA for lens pairs separated 
by i? ± AR/2 as 

N = TV rJg V R AR {Lu + Wgg{R)) as dA (43) 


The galaxy number density in our simulation is Ug = 4.56 x 10 ^/i^Mpc which roughly corresponds to an effective 
absolute magnitude cut of Mr < —20. For this absolute magnitude cut, the SDSS probes an effective survey volume of 


approximately V = 1.2 fsky x 10°/i 


I. 4 or t nis absolute magnitu d 
^Mpc® iTeemark et ahlEoGil . 


where faky is the fraction of the sky covered. For the 


SDSS, fsky will eventually be close to 1/4. The number density of source galaxies used for leasing from isheldon et alJ (l2004h 
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Figure 15. Displayed are logK (top left); I 7 I (top right); 71 (bottom left); 72 (bottom right) The galaxy separation is R12 = 0.8 h ^ 
Mpc and the box size is 4.2 h~^ Mpc on a side. 


is S200/deg^. The area element dA in deg^ is dA=(180 dxj'KD{zLW where dx is the pixel size in physical units that we 
will use for our shear map and D{zl) is the angular diameter distance to redshift zl- We will choose zl = 0.1, the median 
redshift for the SDSS spectroscopic survey, and 7000 deg^ for the estimated sky coverage of the completed SDSS, which gives 
fsky = 0.17. With these numbers our expression for N becomes 

N = 2.0 X leThh Mpc~^ R AR {Lu + Wgg{R)) dx^ (44) 


With this for mula we are ab l e to a dd noise to our shear maps that should correspond to what we will see if we do an 
analysis similar to Isheldon et al ] (I 2 OO 4 I) with SDSS data in the near future to measure the mass density around pairs. We 
take our predicted shear map shown in Figure ESI which corresponds to R12 = 0.8 h ^ Mpc, and add Gaussian noise with 
a = 0.2/FN to each component of the shear. The resulting noisy shear map is shown in Figure El The basic tangential 
pattern can still be seen behind the noise. 

We can now process both the noisy s hear map and the noiseles s one to reconstruct k with a KS mass reconstruction 
algorithm. We use the “direct method” of iLombardi fc BertinI lllOOlil . This inversion routine reconstructs the mass map by 
solving for the Fourier modes on a chosen grid that minimize the functional 




d^ x\\V K,{x) — u(i;)|| 


(45) 


where 


u{x) 


71.1 + 72,1 

72.1 — 71,2 


(46) 


The Euler-Lagrange equations of this functional are just Equations El that we saw related the measured shear to k. This 
routine is linear and can be done quickly with East Eourier Transforms. Unlike the earliest KS algorithms this one does not 
have edge effect biases other than slightly increased noise in the bounding bins. The process involves choosing a grid scale, 
averaging the shear in the grid bins, calculating u{x) by hnite differencing the shear values, and then running Lombardi 
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Figure 16. Displayed are 'yx (top left); •yn (top right); |7t|/|7| (bottom left); |7i{|/i7! (bottom right); The galaxy separation is 
Ri 2 = 0.8 h~^ Mpc and the box size is 4.2 h~^ Mpc on a side. 


and Bertin’s routine (which is publicly available) on the array u. The resulting reconstruction will be degraded in resolution 
for two reasons. First, one will usually want to choose a grid scale that gives reasonable signal-to-noise per pixel. This is 
purely a pixelization effect and structures smaller than the pixel scale will be lost. Second, the reconstruction algorithm which 
differences neighbors introduces an effective smoothing between neighboring pixels. 

We display the reconstruction in Figure IlDl The top right is the original k binned down to the same grid resolution of the 
reconstruction. The top right shows the k recovered from the noisy shear map via the reconstruction algorithm. The bottom 
left shows the residual, which is consistent with zero. The bottom right shows a slice through both maps that intersects both 
peaks around the galaxies. Since the reconstruction cannot recover the constant mass sheet we have subtracted off the mean 
from each map before plotting. 


11 CONCLUSION 

We have calculated the weak lensing shear around pairs of galaxies and related it directly to the 3D galaxy-galaxy-mass 
three-point correlation function. We have shown that the 2D mass profile can be reconstructed from the shear map around 
galaxy pairs and we have discussed how to interpret this in terms of the correlation functions. With our N-body simulation 
and HOD bias model we have calculated the expected signal and have estimated the amount of noise to be expected in a model 
survey similar to the Sloan Digital Sky Survey. Although this prediction is only meant to be a first attempt at estimating the 
effect and does not exhaust the possibilities we predict that this signal may be measurable with current SDSS lensing data. 
The anisotropic dark matter profile about two galaxies extends far away from their center and will produce an anisotropic 
tangential shear profile. Measurements confirming this will be difficult to reconcile with alternative theories of gravity that 
need to produce weak lensing shear without any extended anisotropic dark matter halos. On the other hand we have seen 
that they are a natural prediction of CDM. The three-point function better specifies how mass traces galaxies and so provides 
more detailed information over that contained in the isotropic galaxy-mass two-point correlation function. This exact form of 
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Figure 17. The tangential shear in an annulus of radius R = 2.0 ± 0.1 h~^ Mpc, centered between lens galaxies divided by the mean 
value in the annulus as a function of 9. A quadrupole pattern is evident and we over-plot a pure quadrupole /^ToiP) = (1 + ecos20) 
with and ellipticity parameter e = 0.07. 


the measured galaxy-galaxy-mass correlation function will provide additional constraints on galaxy bias models that govern 
how galaxies trace the underlying dark matter. 
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Figure 19. The top left shows the original k map binned down to the same resolution of the reconstructed one. The top right shows 
the reconstruction obtained by inverting the noisy shear map. The bottom left is the residual from subtracting the top two maps. One 
can see that the residual is consistent with zero but the noise is also correlated over the scale of just a couple of pixels. The bottom right 
plot shows a slice through the two peaks for both densities. One can see that the reconstruction works very well. The means have been 
subtracted off both maps before plotting since this constant mass sheet is unmeasurable. 


Huchra, J. P., Vogeley, M. S., & Geller, M. J. 1999, VizieR Online Data Catalog, 212, 10287 
Hudson, M. J., Gwyn, S. D. J., Dahle, H., & Kaiser, N. 1998, ApJ, 503, 531 
Irgens, R. J., Lilje, P. B., Dahle, H., & Maddox, S. J. 2002, ApJ, 579, 227 

Jenkins, A., Frenk, C. S., Pearce, F. R., Thomas, P. A., Colberg, J. M., White, S. D. M., Couchman, H. M. P., Peacock, 
J. A., Efstathiou, G., & Nelson, A. H. 1998, ApJ, 499, 20 

Jenkins, A., Frenk, G. S., White, S. D. M., Colberg, J. M., Cole, S., Evrard, A. E., Couchman, H. M. P., & Yoshida, N. 2001, 
MNRAS, 321, 372 

Jimenez, R., Thejll, P., Jorgensen, U. G., MacDonald, J., & Pagel, B. 1996, MNRAS, 282, 926 

Joffre, M., Fischer, P., Frieman, J., McKay, T., Mohr, J. J., Nichol, R. C., Johnston, D., Sheldon, E., & Bernstein, G. 2000, 
ApJ, 534, L131 

Johnston, D., Sheldon, E., & Frieman, J. 2004, In preparation 
Kaiser, N. & Squires, G. 1993, ApJ, 404, 441 

Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 1999, MNRAS, 303, 188 
Kauffmann, G. & White, S. D. M. 1993, MNRAS, 261, 921 

Kneib, J.-P., Ellis, R. S., Smail, I., Couch, W. J., & Sharpies, R. M. 1996, ApJ, 471, 643 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

Lombardi, M. & Berlin, G. 1999, A&A, 348, 38 

Luppino, G. A. & Kaiser, N. 1997, ApJ, 475, 20 

Ma, G. & Fry, J. N. 2000a, ApJ, 543, 503 

—. 2000b, ApJ, 531, L87 

McClelland, J. & Silk, J. 1977, ApJ, 217, 331 







26 D. Johnston 


McKay, T. A. et al. 2001 

Milgrom, M. 1983, ApJ, 270, 365 

Natarajan, P. & Refregier, A. 2000, ApJ, 538, LI 13 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Neyman, J. & Scott, E. L. 1952, ApJ, 116, 144 

Peacock, J. A., Cole, S., Norberg, P., Baugh, C. M., Bland-Hawthorn, J., Bridges, T., Cannon, R. D., Colless, M., Collins, 
C., Couch, W., Dalton, G., Deeley, K., De Propris, R., Driver, S. P., Efstathiou, G., Ellis, R. S., Frenk, C. S., Glazebrook, 
K., Jackson, G., Lahav, O., Lewis, L, Lumsden, S., Maddox, S., Percival, W. J., Peterson, B. A., Price, L, Sutherland, W., 
& Taylor, K. 2001, Nature, 410, 169 
Peacock, J. A. & Smith, R. E. 2000, MNRAS, 318, 1144 

Peebles, P. 1980, The Large Seale Structure of the Universe (Princeton University Press) 

Peebles, P. J. E. 1974a, A&A, 32, 197 
—. 1974b, A&A, 32, 197 

Peebles, P. J. E. & Groth, E. J. 1975, ApJ, 196, 1 
Pen, U. 1998, ApJ, 504, 601 

Perlmutter, S., Aldering, G., Goldhaber, G., Knop, R. A., Nugent, P., Castro, P. G., Deustua, S., Fabbro, S., Goobar, A., 
Groom, D. E., Hook, 1. M., Kim, A. G., Kim, M. Y., Lee, J. G., Nunes, N. J., Pain, R., Pennypacker, C. R., Quimby, R., 
Lidman, C., Ellis, R. S., Irwin, M., McMahon, R. G., Ruiz-Lapuente, P., Walton, N., Schaefer, B., Boyle, B. J., Filippenko, 
A. V., Matheson, T., Fruchter, A. S., Panagia, N., Newberg, H. J. M., Couch, W. J., & The Supernova Cosmology Project. 
1999, ApJ, 517, 565 

Press, W. H. & Schechter, P. 1974, ApJ, 187, 425 

Riess, A. G., Filippenko, A. V., Ghallis, P., Glocchiatti, A., Diercks, A., Garnavich, P. M., Gilliland, R. L., Hogan, G. J., Jha, 
S., Kirshner, R. P., Leibundgut, B., Phillips, M. M., Reiss, D., Schmidt, B. P., Schommer, R. A., Smith, R. C., Spyromilio, 
J., Stubbs, C., Suntzeff, N. B., & Tonry, J. 1998, AJ, 116, 1009 
Saunders, W., Rowan-Robinson, M., & Lawrence, A. 1992, MNRAS, 258, 134 
Schneider, P. & Bartelmann, M. 1997, MNRAS, 286, 696 

Scoccimarro, R., Feldman, H. A., Fry, J. N., & Frieman, J. A. 2001a, ApJ, 546, 652 

Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001b, ApJ, 546, 20 

Scranton, R. 2003, MNRAS, 339, 410 

Seitz, S. & Schneider, P. 2001, A&A, 374, 740 

Seljak, U. 2000, MNRAS, 318, 203 

Sheldon, E., Johnston, D., Frieman, J., Scranton, R., McKay, T., Gonnolly, A., & SDSS collaboration. 2004, In preparation 
Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119 

Smith, D. R., Bernstein, G. M., Fischer, P., & Jarvis, M. 2001, ApJ, 551, 643 

Spergel, D. N., Verde, L., Peiris, H. V., Komatsu, E., Nolta, M. R., Bennett, C. L., Halpern, M., Hinshaw, G., Jarosik, N., 
Kogut, A., Limon, M., Meyer, S. S., Page, L., Tucker, G. S., Weiland, J. L., Wollack, E., & Wright, E. L. 2003, ApJS, 148, 
175 

Takada, M. & Jain, B. 2003a, MNRAS, 340, 580 

—. 2003b, MNRAS, 344, 857 

Tegmark, M. et al. 2003 

Totsuji, H. & Kihara, T. 1969, PASJ, 21, 221 

Tyson, J. A. & Fischer, P. 1995, ApJ, 446, L55+ 

Tyson, J. A., Valdes, F., Jarvis, J. F., & Mills, A. P. 1984, ApJ, 281, L59 
White, M. 2002, ApJS, 143, 241 
White, M. J. 2003 

White, S. D. M. & Frenk, C. S. 1991, ApJ, 379, 52 

White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 

Wilson, G., Kaiser, N., Luppino, G. A., & Gowie, L. L. 2001, ApJ, 555, 572 

York, D. G. & SDSS collaboration. 2000, AJ, 120, 1579 

Zehavi, 1. et al. 2003, ApJ, 608 ,16 



